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TTTSAMME ^ASSUNG 
Computertomographieverfahren mit kegelformigen Strahlenbundeln 



Die Erfindung betrifft ein Computertomographieverfahren, bei dem ein Rebinning zu 
parallelen Strahlenfachern erfolgt, deren Strahlen eine die Rotationsachse enthaltene 
Ebene in aquidistanten DurchstoBpunkten durchsetzen. Die Rekonstruktion erfolgt fur 
Strahlen, die auf diesen Ebenen senkrecht stehen. Durch diesen Ubergang von einer 
Kegelstrahl-Geometrie zu einer Parallelstrahl-Geometrie ist eine sehr schnelle 
Rekonstruktion moglich, die insbesondere in der CT-Fluoroskopie Anwendung finden 
kann. 



Fig. 6 




FIG. 6 



PHDE020092 



BESCHREffiUNG 
Computertomographieverfahren 

Die Erfindung betrifft ein Computertomographieverfahren, bei dem ein Untersuchungs- 
bereich mit einem kegelformigen Strahlenbundel durchstrahlt wird und bei dem die 
Strahlenquelle, die das Strahlenbundel erzeugt, relativ zum Untersuchungsbereich auf 
einer Kreisbahn urn eine Rotationsachse rotiert. Aufierdem betrifft die Erfindungen einea 
Computertomographen zur DurchfUhrung des Verfahrens sowie ein Computerprogramm 
zur Steuerung des Computertomographen. 

Ein Verfahren dieser Art und ein entsprechender Computertomograph ist aus der 
US PS 6,285,733 (PHD 98-123) bekannt Die dabei von einer zweidimensionalen 
Detektoreinheit akquirierten Messwerte, die von der Intensitat in dem Strahlenbundel 
jenseits des Untersuchungsbereiches abhangen, werden zunachst einem Rebinning 
unterzogen. Dabei werden Gruppen von Messwerten gebildet, die zu Strahlenfachern 
gehoren, die sich in zueinander und zur Rotationsachse parallelen, aquidistanten 
Facherebenen befinden. Die Strahlenfacher setzen sich nach dem Rebinning aus Strahlen 
zusammen, die eine die Rotationsachse enthaltende, zu den Facherebenen der jeweiligen 
Gruppe senkrechte Ebene in DurchstoBpunkten durchsetzen, die sich auf zur Rotations- 
achse senkrechten und zueinander parallelen, aquidistanten Verbindungsgeraden 
befinden. 

Die durch das Rebinning entstandenen Messdaten werden anschlieflend einer 
eindimensionalen Hochpassfilterung sowie einer Rttckprojektion zwecks Erzeugung 
wenigstens eines CT-Bildes unterzogen. Dabei werden zur Rekonstruktion eines Voxels 
im Untersuchungsbereich von jeder Gruppe die Messdaten von Strahlen herangezogen, 
die das betreffende Voxel aus unterschiedlichen Richtungen durchsetzt haben. Wenn sich 
das Voxel nicht in der durch die Kreisbahn defmierten Zentralebene befindet, verlaufen 
die Strahlen durch das Voxel unter einem Winkel zu der Zentralebene, und jeder Strahl 
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durchsetzt die zu seiner Gruppe gehorende Ebene unter einem von 90° verschiedenen 
Winkel. Die Rekonstruktion unter Berucksichtigung dieser Kegelstrahl-Geometrie ist 
relativ kompliziert und erfordert dementsprechend relativ viel Rechenzeit. 

Bei verschiedenen Anwendungen sind aber verkurzte Rechenzeiten erwunscht. 
Beispielsweise bei der CT-gefiihrten Biopsie (CT = Computertomographie), wobei eine 
Biopsienadel in ein Untersuchungsobjekt eingefuhrt wird und der Vorschub der 
Biopsienadel fortlaufend anhand einer Folge von dreidimensionalen CT-Bildern 
kontrolliert wird. Bei dieser sogenannten CT-Fluoroskopie steht fiir die Rekonstruktion 
eines CT-Bildes nur sehr wenig Zeit zur Verfiigung, insbesondere wenn dabei eine 
sogenannte „Sliding-Window"-Technik benutzt wird, bei der ein Update eines CT-Bildes 
(unter Heranziehung von neu akquirierten CT-Daten und von CT-Daten, die bereits bei 
der Rekonstruktion des vorigen CT-Bildes benutzt wurden) innerhalb eines Zeitraums 
erfolgt, der wesentlich kurzer ist als der fur die Acquisition der Messdaten eines 
vollstandigen CT-Bildes benotigte Zeitraum. 

Aufgabe der vorliegenden Erfindung ist es daher, ein Verfahren der eingangs genannten 
Art so auszugestalten, dass die fur die Rekonstruktion eines CT-Bildes benotigte Zeit 
verringert wird> Diese Aufgabe wird erfindungsgemaB gelost durch ein 
Computertomographie- Verfahren mit den Schritten: 

a) Erzeugen eines kegelfbrmigen, einen Untersuchungsbereich bzw. ein darin 
befindliches Objekt durchsetzenden Strahlenbundels mit einer Strahlenquelle, 

b) Erzeugung einer kreisformigen, eine Rotation urn eine Rotationsachse umfassenden 
Relativbewegung zwischen der Strahlenquelle einerseits und dem 

25| Untersuchungsbereich bzw. dem Objekt andrerseits, 

c) Akquisition von Messwerten, die von der Intensitat in dem Strahlenbiindei jenseits 
des Untersuchungsbereiches abhangen, mit einer Detektoreinheit wahrend der 
Relativbewegung, 

d) Rebinning der Messwerte zu einer Anzahl von Gruppen, wobei jede Gruppe die 
Messwerte von Strahlenfkchern umfasst, die sich in zueinander und zur 
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Rotationsachse parallelen, aquidistanten Facherebenen befinden und die sich aus 
Strahlen zusammensetzen, die eine die Rotationsachse enthaltende, zu den 
Facherebenen dieser Gruppe senkrechte Ebene in DurchstoBpunkten durchsetzen, 
die sich auf zur Rotationsachse senkrechten und zueinander paralleled aquidistanten 
; > Verbindungsgeraden befinden, 

e) Rekonstruktion der raumlichen Verteilung der Schwachung der Rontgenstrahlung 
aus den durch das Rebinning der Messwerte gebildeten Messdaten fur Strahlen, die 
auf den Ebenen der Gruppen senkrecht stehen und die Durchstofipunkte durchsetzen 
zur Erzeugung wenigstens eines CT-Bildes. 

10 

Wahrend bei dem bekannten Verfahren die Rekonstruktion unter Berucksichtigung der 
bei der Akquisition der Messwerte gegebenen Kegelstrahl-Geometrie erfolgt, wird der 
Rekonstruktion bei der Erfindung eine Parallelstrahl-Geometrie zugrunde gelegt, bei der 
alle Strahlen senkrecht zu der Ebene der zugehorigen Gruppe und damit parallel zu der 
IS durch die Kreisbahn definierten Zentralebene verlaufen. Der Rechenaufwand fur eine 
solche Parallelstrahl-Geometrie ist wesentlich geringer als fur eine Kegelstrahl- 
Geometrie. Ein weiterer Vorteil ist, dass bei der Erfindung ein zylinderfbrmiger 
Abschnitt des Untersuchungsbereichs rekonstruiert wird. 

20 Da bei der Erfindung der Kegelwinkel der Strahlen (das ist der Winkel, den die Strahlen 
mit einer zur Rotationsachse senkrechten Ebene einschlieBen) de facto ignoriert wird, 
lasst sich der Untersuchungsbereich aufierhalb der Zentralebene nicht exakt rekon- 
struieren, was - wenn auch aus anderen Griinden - auch fur das bekannte Verfahren gilt. 
Der Verlust an Bildqualitat ist aber gering, wenn die Detektoreinheit nur wenige, 

2f Detektorzeilen umfasst bzw. wenn der maximale Kegelwinkel gering ist. Insbesondere ist 
die Bildqualitat besser als bei einem Verfahren, bei dem von vornherein die Messwerte 
so behandelt werden als wlirden sie in parallelen Ebenen verlaufen, in denen sich jeweils 
eine Detektorzeile befindet 

30 
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Die Anspriiche 2 und 3 definieren bewahrte Verfahren zur Rekonstruktion eines CT- 
Bildes bei einer ParallelstraW-Geometrie. Diese Verfahren erlauben an sich eine exakte 
Rekonstruktion des Untersuchungsbereiches, doch ist die Rekonstruktion nur innerhalb 
der Zentralebene exakt, weil die der Rekonstruktion zugrunde liegende Parallelstrahl- 
Geometrie bei der Akquisition der Messwerte nur in der Zentralebene gegeben ist. 

Anspruch 4 definiert eine bevorzugte Anwendung des erfindungsgemafien Verfahrens; 
die CT-Fluorosko P ie erfordert eine schnelle Rekonstruktion der CT-Bilder, insbesondere 
wenn die zeitlichen Abstande zwischen zwei Updates eines CT-Bildes kurzer sind als der 
id Zeitraum, der fur die vollstandige Akquisition der Messwerte fur ein CT-Bild benotigt 
wird. 

Anspruch 5 beschreibt einen Computertomographen zur Durchfuhrung des erfindungs- 
gemaBen Verfahrens und Anspruch 6 ein Computertomogramm zur Steuerung eines 
15| solchen Computertomographen. 

Die Erfindung wird nachstehend anhand der Zeichnungen naher erlautert. Es zeigen: 
Fig. 1 einen Computertomographen, mit dem die Erfindung ausftihrbar ist, 
Fig. 2 ein Ablaufdiagramm fur das erfindungsgemaBe Verfahren, 
Fig. 3 ein in einer Strahlenquellenposition erzeugtes kegelfdrmiges Strahlenbundel, 
Fig. 4 die durch das Rebinning gebildeten Strahlenfacher in parallelen Ebenen, 
Fig. 5 einen Querschnitt durch diese Strahlenfacher, 
Fig. 6 die Kegelstrahl- und die Parallelstrahl-Geometrie. 



1 
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Der in Fig. 1 dargestellte Computertomograph umfasst eine Gantry 1, die urn eine 
parallel zur z-Richtung des in Fig. 1 dargestellten Koordinatensystems verlaufende 
Rotationsachse 14 rotieren kann. Dazu wird die Gantry von einem Motor 2 mit einer 
vorzugsweise konstanten aber einstellbaren Winkelgeschwindigkeit angetrieben. An der 
Gantry ist eine Strahlenquelle S, beispielsweise ein Rdntgenstrahler, befestigt. Dieser ist 
mit einer Kollimatoranordnung 3 versehen, die aus der von der Strahlenquelle S 
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erzeugten Strahlung ein kegelformiges Strahlenbundel 4 ausblendet, d. h. ein Strahlen- 
bundel, das sowohl in z-Richtung als auch in einer dazu senkrechten Richtung (d. h. in 
einer zur Rotationsachse senkrechten Ebene) eine von Null verschiedene endliche 
Ausdehnung hat. 

Das Strahlenbundel durchdringt einen Untersuchungsbereich 13, in dem sich ein Unter- 
suchungsobjekt, z. B. ein Patient auf einem Patientenlagerungstisch (beides nicht naher 
dargestellt), befinden kann. Der Untersuchungsbereich 13 hat die Form eines Zylinders. 
Nach dem Durchsetzen des Untersuchungsbereichs 13 trifift das R6ntgenstrahlenbundel 4 
auf eine an der Gantry 1 befestigte zweidimensionale Detektoreinheit 16. Diese umfasst 
eine Anzahl von in z-Richtung nebeneinander angeordneten Detektorzeilen mit jeweils 
einer Vielzahl von Detektorelementen. Die Detektorzeilen befinden sich in zur 
Rotationsachse senkrechten Ebene auf einem Kreisbogen urn die Strahlenquelle S; sie 
konnen aber auch einen Kreisbogen urn die Rotationsachse 14 beschreiben oder 
geradlinig sein. Jedes Detektorelement liefert in jeder Strahlenquellenposition einen 
Messwert for einen Strahl aus dem Strahlenbundel 4. 

Der mit cwbezeichnete Offhungswinkel des Strahlenbundels (als Offhungswinkel ist der 
Winkel definiert, den ein Strahl, der in einer zur Rotationsachse 14 senkrechten Ebene ... 
am Rande des Strahlenbundels 4 liegt, mit einer durch die Strahlenquelle S und die 
Rotationsachse 14 definierten Zentralstrahlebene einschlieBt) bestimmt den Durchmesser 
des Zylinders, innerhalb dessen sich das zu untersuchende Objekt bei der Akquisition der 
Messwertebefindet. Der Untersuchungsbereich 13 bzw. das Untersuchungsobjekt oder 
der Patientenlagerungstisch kann mittels eines Motors 5 parallel zur Rotationsachse 14 
verschoben werden. Wenn die beiden Motoren 2 und 5 gleichzeitig laufen, ergibt sich 
eine helixformige Bewegung der Strahlenquelle S und der Detektoreinheit 16 urn den 
Untersuchungsbereich 13. Wenn hingegen der Motor 5 fiir den Vorschub in z-Richtung 
stillsteht und nur der Motor 2 die Gantry einzeln rotieren lasst, ergibt sich. eine 
kreisfbrmige Abtastwegung der Strahlenquelle S und der Detektoreinheit 16 relativ zum 
Untersuchungsbereich 13. 
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Die von der Detektoreinheit 16 akquirierten Messdaten werden von einem Bildver- 
arbeitungsrechner 10 zugefuhrt, der daraus die Absorptionsverteilung in einem Teil des 
Untersuchungsbereichs 13 rekonstruiert und z. B. auf einem Monitor wiedergibt Die 
beiden Motoren 2 und 5, der Bildverarbeitungsrechner 10, die Strahlenquelle S und der 
Transfer der Messwerte von der Detektoreinheit 16 zum Bildverarbeitungsrechner 10 
werden von einer geeigneten Kontrolleinheit 7 gesteuert. 

Fig. 2 zeigt den Ablauf eines Akqusitions - und Rekonstruktionsverfahrens, das mit dem 
ComputertomographennachFig. 1 durchgefuhrt -werden kann. 

Nach der Initialisierung im Block 101 rotiert die Gantry mit einer konstanten 
Geschwindigkeit, wobei sich eine Umlaufzeit von 1 sec oder weniger ergibt. Die 
Strahlenquelle S emittiert ein kegetf6rmiges Strahlenbundel, das den Untersuchungs- 
bereich durchsetzt, und die von den Detektorelementen der Detektoreinheit 16 
akquirierten Messwerte werden im Bildverarbeitungsrechner 10 zwischengespeichert 
und weiterverarbeitet. 

Fig. 3 zeigt die kreisformige Bahn 17, auf der sich die Strahlenquelle S und die 
Detektoreinheit 16 urn die Rotationsachse 14 bewegen. Die Bahn 17 liegt in einer zur 
Rotationsachse senkrechten Ebene, die im folgenden als Zentralebene bezeichnet wird. 
Es ist ein Strahlenbundel 4 dargestellt, das in einer bestimmte Strahlenquellenposition So 
von der Strahlenquelle emittiert wird. Dieses kegelfbrmige Strahlenbundel kann in eine 
Vielzahl von ebenen Strahlenfachern zerlegt werden, die sich - wie die in Fig. 3 
dargestellten Strahlenfacher 401, 402, 403 - in zur Rotationsachse 14 parallelenJEbenen 
befinden. Die Strahlenfacher gehen von derselben Strahlenquellenposition aus und 
werden von je einer zur Rotationsachse 14 parallelen Spalte von Detektorelementen auf 
der Detektoreinheit 16 erfasst. 



3C 



Fig. 3 deutet an, dass auch in anderen Positionen der Strahlenquelle, z. B. S. 2> S.,, Si 
oder S 2 , das emittierte kegelfbrmige Strahlenbundel gemessen wird. Diese Strahlen- 




PHDE020092 

-7- 

quellenpositionen bzw. die dort emittierten Strahlenbundel konnen durch einen 
Parameter P charakterisiert werden, der dem Winkel entspricht, den das Lot von der 
Strahlenquellenposition auf die Rotationsachse 14 mit einer Referenzgeraden in der 
Zentralebene einschlieBt (0 kann entsprechend der Anzahl der Umdrehungen der 
Strahlenquelle urn die Rotationsachse groBer sein als 2n). Die Lage eines jeden 
Strahlenfachers in einem Strahlenbundel kann durch einen Parameter s gekennzeichnet 
werden, der die Lage der von dem Strahlenfacher getroffenen Spalte von Oetektor- 
elementen innerhalb der Detektoreinheit 16 beschreibt. Jeder Strahl innerhalb eines 
solchen Strahlenfdchers kann wiederum durch den Parameter t charakterisiert werden, 
der die Lage des von dem betreffenden Strahl getroffenen Detektorelements innerhalb 
der Spalte von Detektorelementen beschreibt bzw. den Abstand dieses Dektorelements 
von der Zentralebene. 



Somit bilden die akquirierten Messwerte einen dreidimensionalen Datensatz M 0 (P, s, t), 
s wobei jeder Messwert einem Gitterpunkte eines regelmafligen kartesischen Gitters in 
einem dreidimensionalen (P, s, t)-Parameterraum entspricht. Die Akquisition der 
Messwerte stellt daher eine Abtastung der sogenannten Objektfonktion (in diesem Fall 
des Linienintegrals der Schwachung der Strahlung) an einer Vielzahl von gleichmaBig im 
(P, s, t)- Parameterraum verteilten Punkten dar. 
20 

Die Akquisition der Messwerte im Schritt 102 und die Verarbeitung dieser Messwerte in 
den darauf folgenden Schritten 103 ff erfolgen zeitlich parallel zueinander, so dass die 
vorhandenen Messwerte bereits weiterverarbeitet werden, wahrend weitere Messwerte 
akquiriert werden. 

25 

Im Schritt 103 werden die Messwerte mit dem Cosinus des Winkels multipliziert, den 
der Strahl, entlang dessen der Messwert erfasst wurde, mit einer zur Rotationsachse 
senkrechten Ebene einschlieBt. Wenn die Abmessungen der Detektoreinheit in z- 
Richtung gering sind, kann dieser Schritt aber entfallen, weil der Winkel dann so klein 
3C ist, dass der Cosinus des Winkels stets nahezu 1 ist. 
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Der im Schritt 102 akquirierte und im Schritt 103 gegebenenfalls modifizierte Datensatz 
Mo(P, s, t) ist fur die weitere Verarbeitung noch nicht optimal. Deshalb erfolgt in den 
Schritten 104 und 105 ein sogenanntes Rebinning der Messwerte. Die Daten werden 
dabei so umsortiert und uminterpoliert, als waren sie mit einer anderen Strahlenquelle 
(einer kreisbogenfdrmigen Strahlenquelle, die zueinander parallele Strahlenfacher 
emittiert) und mit einem anderen Detektor (einem ebenen, rechteckigen und die 
Rotationsachse enthaltenden Detektor) gemessen worden. In Schritt 104 werden zu 
diesem Zweck zunachst die in zueinander und zur Rotationsachse 14 parallelen Ebenen 
liegenden Strahlenfacher von verschiedenen Strahlenquellenpositionen zu jeweils einer 
1(| Gruppe zusammengefasst. 

Fig. 4 zeigt eine auf diese Weise gebildete Gruppe von Strahlenfachern. Von jeder 
Strahlenquellenposition S* ... S 0 ... S 2 gehort jeweils ein Strahlenfacher zu dieser 
Gruppe. Die zu einer Gruppe gehorenden Strahlenfacher erfullen die Bedingung 

1^ <j)=a+p _ (1) 

Dabei ist <j> die Projektionsrichtung, in der eine Gruppe von Strahlenfachern den 
Untersuchungsbereich durchsetzt. a ist der Winkel, den der betreffende Strahlenfacher in 
dem ursprttnglichen Strahlenbundel (vgl. Fig. 3) mit einer durch die Rotationsachse 14 
und die Strahlenquellenposition (z. B. S 2 , die ihrerseits durch den Winkel p definiert ist) 
2(| defmierten Ebene einschlieBt. Derartige Gruppen von Strahlenfachern werden fur 
verschiedene Projektionsrichtungen $ gebildet, die sich jeweils um ein bestimmtes 
Projektionswinkel-Inkrement d<|) voneinander unterscheiden. Wenn die Strahlenfacher 
eines Strahlenbiindels die Gleichung 1 nicht exakt erfullen, muss ein entsprechender 
Strahlenfacher durch Interpolation aus den Strahlen benachbarter Strahlenfacher des 
Strahlenbiindels 4 (Fig. 3) ermittelt werden. 



25 
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Die Strahlenfacher einer Gruppe - darunter die in Fig. 4 dargestellten Strahlenfacher 41 1 
... 415 - definieren ein Strahlenbundel 410, das eine zeltartige Form hat und sich aus 
Strahlenfachern zusammensetzt, die in zueinander und zur Rotationsachse parallelen 
Ebenen liegen. In Fig. 4 ist aufJerdem der Schnittbereich 420 dargestellt, der sich ergibt, 
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wenn das Strahlenbundel 410 von einer die Rotationsachse 14 enthaltenden und zu den 
Ebenen der Strahlenfacher 41 1 ... 415 senkrechten Ebene geschnitten wird. Die oberen 
und unteren Rander sind kissenformig nach auBen gewolbt, weil die Strahlenquellen- 
positionen in der Mitte weiter von der Schnittebene entfernt liegen als die Strahlen- 
quellenpositionen am Rande. 

Fig. 5 stellt den Schnittbereich 420 fan Detail dar. Durch runde Punkte sind die Stellen 
markiert, an denen die Strahlen der Strahlenfacher 411 ... 415 den Schnittbereich 
durchstoBen. Diese Punkte liegen bei den auBeren Strahlenfachern 411 bzw. 415 enger 
beisammen als bei den inneren Strahlenfachern 412 ... 414, und in einem der inneren 
Strahlenfacher (z. B. 413) nimmt der Abstand dieser Punkte voneinander von der Mitte 
nach auBen hin ab. In dem Schnittbereich 420 ist ein rechteckiger Bereich 160 enthalten, 
dessen oberer und unterer Rand 161 bzw. 162 durch die Abmessungen der beiden 
auBersten Strahlenfacher 41 1 bzw. 415 im Schnittbereich 420 gegeben sind und der im 
Folgenden auch als „virtueller Detektor" bezeichne wird. 

Im Schritt 105 Werden die Strahlen der Strahlenfacher 41 1 ... 415 einer Gruppe durch 
Interpolation neu berechnet, und zwar so, dass sie den virtuellen Detektor 160 in 
DurchstoBpunkten durchsetzen, die sich auf zur Rotationsachse 14 senkrechten und 
. zueinander parallelen, aquidistanten Verbindungsgeraden befinden. Auf jeder Ver- 
bindungsgeraden ist der Abstand der DurchstoBpunkte voneinander konstant. Die aus 
diesem Rebinningschritt resultierenden DurchstoBpunkte sind in Fig. 5 durch ein „+" 
markiert, und eine der erwahnten, zu den Seiten 161 und 162 parallelen Verbindungs- 
geraden ist in Fig. 5 durch die gestrichelte Gerade G dargestellt. 

Durch die beiden Rebinningschritte 104 und 105 werden also die Messwerte Mo(B, s, t), 
die durch ein regelmaBiges Gitter im (p, s, t)-Parameterraum definiert sind, uberfuhrt in 
Messwerte Mi(<|>, u, v), die durch ein regelmaBiges Gitter in einem dreidimensionalen 
u, v)-Parameterraum definiert sind. Die Parameter u bzw. v stellen die Koordinaten der 
DurchstoBpunkte senkrecht bzw. parallel zur Rotationsachse in dem virtuellen Detektor 
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160 dar. Die durch das Rebinning entstandenen Werte Mi(*. u. v) werden im Folgenden 
als Messdaten bezeichnet. um sie begrifflich von den aus der Messung hervorgehenden 
MesswertenMo(P, s, t) zu unterscheiden. 

Da die Strahlen, die den Schnittbereich 420 oberhalb des oberen Randes 161 bzw. 
unterhalb des unteren Randes 162 durchsetzen, im weiteren Verfahren nicht mehr 
benutzt werden, ist es von Vorteil die Kollimatoranordnung 3 (Fig. 1) so zu gestalten, 
dass das kegelformige Strahlenbundel diese Strahlen nicht enthalt. Anstelle von 
geradlinigen senkrecht zur Rotationsachse verlaufenden Kanten miisste die Kollimator- 
anordnung 3 zu diesem Zweck nach innen gewolbte Kanten aufweisen. Dadurch wurde 
die Strahlenbelastung fur den Patienten verringert. 

Insoweit als bisher anhand von Fig. 2 erlautert, ist die Erfindung aus der US PS 6285733 
(PHD 98-123) bekannt. ErfindungsgemaB erfolgt die weitere Verarbeitung der Mess- 
daten aber nicht auf Basis einer Kegelstrahl-Geometrie, sondern auf Basis einer ParaUel- ■ 
strahl-Geometrie. Dies ist durch den gestrichelten Kasten 106 symbolisch angedeutet 
und in Fig. 6 veranschaulicht. Fig. 6 zeigt mit ausgezogenen Linien die divergierenden, 
von einer Strahlenquellenposition (z. B. So) ausgehenden Strahlen eines Strahlenfachers 
(z. B. 413). Diese Strahlen durchsetzen den Virtuellen Detektor 160 in aquidistanten 
Durchstofipunkten. Die Erfindung verwendet stattdessen die gestrichelt dargestellten 
Strahlen, die senkrecht auf den virtuellen Detektor 160 auftreffen. Die durch das 
Rebinning entstandenen Messdaten werden den parallelen Strahlen zugeordnet, die 
jeweils denselben DurchstoBpunkt des virtuellen Detektors 160 durchsetzen wie die von 
der Strahlenquellenposition S 0 ausgehenden divergenten Strahlen. 

Dieser Ubergang zu parallelen Strahlen erfordert keinen gesonderten Rechenschritt 
(weshalb der Kasten 106 auch nur gestrichelt dargestellt ist), weil es lediglich notwendig 
ist, den Kegelwinkel (das ist der Winkel, den die Strahlen mit einer Rotationsachse 14 
senkrechtenEbene einschlieBen) zu ignorieren. Uber einen Mehrfachschichtverfahren 
(Multi-Slice-Verfahren), bei dem Daten in parallelen Schichten akquiriert werden, ist 





PHDE020092 



- 11 - 



dann eine Rekonstruktion dadurch m6gtich, dass unter Heranziehung der Messdaten von 
Strahlen die jeweils in derselben, zur Rotationsachse 14 senkrechten Ebene liegen, die 
Schichten voneinander unabhangig rekonstruiert werden. 

5 Mit der Erfindung kann somit ein zylinderformiger Teil des Untersuchungsbereichs 
rekonstruiert werden, dessen H6he dem Abstand der beiden Render 161 und 162 des 
virtuellen Fensters 160 voneinander entspricht. Fur die Hohe h des zylinderformigen 
Teils des Untersuchungsbereiches gilt die Beziehung 



Dabei ist r der Radius der Kreisbahn 17, Yma* ist der maximale Kegelwinkel und a ist der 
Facherwinkel, d. h. der Offhungswinkel des Strahlenbundels in einer zur Rotationsachse 
senkrechten Ebene. 



Ausgehend von den geschilderten geometrischen Verhaltnissen erfolgt in den Schritten 
107 bis 109 eine Rekonstruktion durch eine gefilterte Ruckprojektion. Im Schritt 107 
erfolgt eine eindimensioriale rampenfbrmige Filterung der Messdaten Mi(9, u; v). Dabei 
werden alle Messdaten mit demselben Projektionswinkel 0 und mit demselben Parameter 

20 v einer Filterung unterzogen, bei der der Ubertragungsfaktor rampenfbrmig mit der 
Frequenz ansteigt 

Die auf diese Weise gefilterten Daten F(9, u, v) werden dann einer Ruckprojektion 
unterzogen. Dabei wird im Schritt 108 ein Voxel P(x, y, z) vorgegeben. Im Schritt 109 

21 werden danach von den Messdaten, deren zugehdrige Strahlen das Voxel aus ver- 
schiedenen Projektionsrichtungen $ durchsetzen und in derselben, zur Zentralebene 
parallelen Ebene liegen wie das Voxel, BeitrSge zu dem Schwachungswert des 
betreffenden Voxels ermittelt und summiert. Dies wird fur alle Voxel in der erwahnten 
Ebene wiederholt so dass die Objektfunktion flir diese Ebene bekannt ist. Die Schritte 

3( 108, 109 werden fur alle zur Zentralebene parallelen Ebenen und die in' diesen Ebenen 



10 



h=rtan(y majc )cos(a) 



(1) 
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jeweils befindlichen Voxel und Strahlen ausgeftihrt. Danach kann das auf diese Weise 
rekonstruierte CT-Bild dargestellt werden (Schritt 110). 

Bei einem fluoroskopischen Computertomographieverfahren sind die Schritte 101 bis 
5 110 Teile eines kontinuierlich verlaufenden Akquisitions- und Rekonstruktionsprozesses, 
bei dem fortlaufend neue Messwerte akquiriert und neue CT-Bilder des Untersuchungs- 
bereichs erstellt werden. Der zeitliche Abstand zwischen Updates eines CT-Bildes kann 
dabei kurzer sein als die Zeit, die der Computertomograph benotigt, um Messwerte fiir 
ein vollstandiges CT-Bild zu akquirieren. In diesem Fall erfolgt die Rekonstruktion des 
10 CT-Bildes zu einem Teil aus Messwerten, die neu akquiriert wurden, und zu einem 

anderen Teil aus Messwerten, die bereits bei der Rekonstruktion des vorangehenden CT- 
Bildes benutzt wurden. 

Anstatt durch eine gefilterte Ruckprojektion in den Schritten 107 bis 109 kann die 
15 raumliche Verteilung der Schwachung in den durch die DurchstoBpunkten definierten 
Ebenen auch mit Hilfe einer inversen Fourier-Transformation erfolgen, nachdem zuvor 
die Messdaten durch ein Gridding-Verfahren auf ein kartesisches Gitter in den einzelnen 
Ebenen transformiert wurde. 



20 
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PATENTANSPROCHE 



1. Computertomographie-Verfahren mit den Schritten: 

a) Erzeugen eines kegelfdrmigen, einen Untersuchungsbereich bzw. ein darin 
befindliches Objekt durchsetzenden Strahlenbiindels mit einer Strahlenquelle, 

b) Erzeugung einer kreisformigen, eine Rotation urn eine Rotationsachse umfassenden 
Relativbewegung zwischen der Strahlenquelle einerseits und dem 
Untersuchungsbereich bzw. dem Objekt andrerseits, 

c) Akquisition von Messwerten, die von der Intensitat in dem Strahlenbundel jenseits 
des Untersuchungsbereiches abhangen, mit einer Detektoreinheit wahrend der 
Relativbewegung, 

d) Rebinning der Messwerte zu einer Anzahl von Gruppen, wobei jede Gruppe die 
Messwerte von Strahlenfachern umfasst, die sich in zueinander und zur 
Rotationsachse parallelen, aquidistanten Facherebenen befinden und die sich aus 
Strahlen zusammensetzen, die eine die Rotationsachse enthaltende, zu den 
Facherebenen dieser Gruppe senkrechte Ebene in DurchstoBpunkten durchsetzen, 
die sich auf zur Rotationsachse senkrechten und zueinander parallelen, aquidistanten 
Verbindungsgeraden befinden, 

e) Rekonstruktion der raumlichen Verteilung der Schwachung der R6ntgenstrahlung 
aus den durch das Rebinning der Messwerte gebildeten Messdaten fUr Strahlen, die 
auf den Ebenen der Gruppen senkrecht stehen und die DurchstoBpunkte durchsetzen 
zur Erzeugung wenigstens eines CT-Bildes. 
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2. Computertomographie-Verfahren nach Anspruch 1, wobei sich der 
Rekonstruktionsschritt aus den folgenden Schritten zusammensetzt: 

a) Eindimensionale Filterung der durch das Rebinning erzeugten Messdaten einer jeden 
Gruppe in Richtung der Verbindungsgeraden, 

b) Rttckprojektion der gefilterten Daten von mehreren Gruppen. 
r 

3. Computertomographie-Verfahren nach Anspruch 1, wobei der Rekonstruktionsschritt 
eine inverse Fouriertransformation umfasst. 

1<|> 4. Verwendung des Verfahrens nach 1 bei der CT-Fluoroskopie, wobei unter 
fortlaufender Rekonstruktion von CT-Bilder eine fortlaufende Akquisition von 
Messwerten fur weitere CT-Bilder erfolgt. 

5. Computertomograph zur Durchfuhrung des Verfahrens nach Anspruch 1 mit 
1 $ a) einer Strahlenquelle zum Erzeugen eines kegelformigen, einen Untersuchungsbereich 
bzw. ein darin befindliches Objekt durchsetzenden Strahlenbiindels, 
b) einer Antriebsanordnung zur Erzeugung einer kreisfrirmigen, eine Rotation urn eine 
Rotationsachse umfassenden Relativbewegung zwischen der Strahlenquelle 
einerseits und dem Untersuchungsbereich bzw. dem Objekt andrerseits, 
2(j c) einer Detektoreinheit zur Akquisition von Messwerten, die von der Intensitat in dem 
Strahlenbiindel jenseits des Untersuchungsbereiches abhangen, wahrend der 
Relativbewegung, 

und mit einer Bildverarbeitungseinheit zur Erzeugung mindestens eines CT-Bildes aus 
| den Messwerten mit den Schritten: 
2^ d) Rebinning der Messwerte zu einer Anzahl von Gruppen, wobei jede Gruppe die 

I Messwerte von Strahlenfachern umfasst, die sich in zueinander und zur 

i - 

Rotationsachse parallelen, aquidistanten Facherebenen befinden und die sich aus 

! Strahlen zusammensetzen, die eine die Rotationsachse enthaltende, zu den 

; Facherebenen dieser Gruppe senkrechte Ebene in DurchstoBpunkten durchsetzen, 
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die sich auf zur Rotationsachse senkrechten und zueinander paralleled aquidistanten 
Verbindungsgeraden befinden, 
e) Rekonstruktion der raumlichen Verteilung der Schwachung der Rontgenstrahlung 
aus den durch das Rebinning der Messwerte gebildeten Messdaten for Strahlen, die 
auf den Ebenen der Gruppen senkrecht stehen und die DurchstoBpunkte durchsetzen 
zur Erzeugung wenigstens eines CT-Bildes. 

6. Computerprogramm zur Steuerung eines Computertomographen nach Anspruch 5 
gemaB folgendem Ablauf: 
1C a) Erzeugen eines kegelformigen, einen Untersuchungsbereich bzw. ein darin 

befindliches Objekt durchsetzenden Strahlenbundels mit einer Strahlenquelle, 
b) Erzeugung einer kreisformigen, eine Rotation urn eine Rotationsachse umfassenden 
Relativbewegung zwischen der Strahlenquelle einerseits und dem 
Untersuchungsbereich bzw. dem Objekt andrerseits, 
1 5 c) Akquisition von Messwerten, die von der Intensitat in dem Strahlenbiindel jenseits - - 
des Untersuchungsbereiches abhangen, mit einer Detektoreinheit wahrend der 
Relativbewegung, 

d) Rebinning der Messwerte zu einer Anzahl von Gruppen, wobei jede Gruppe die 
Messwerte von Strahlenfachern umfasst, die sich in zueinander und zur 
20 Rotationsachse parallelen, aquidistanten Facherebenen befinden und die sich aus 
j Strahlen zusammensetzen, die eine die Rotationsachse enthaltende, zu den 
| Facherebenen dieser Gruppe senkrechte Ebene in DurchstoBpunkten durchsetzen, 
| die sich auf zur Rotationsachse senkrechten und zueinander parallelen, aquidistanten 

j Verbindungsgeraden befinden, 

i 

25\ e) Rekonstruktion der raumlichen Verteilung der Schwachung der Rdntgenstrahlung 
i aus den durch das Rebinning der Messwerte gebildeten Messdaten flir Strahlen, die 
auf den Ebenen der Gruppen senkrecht stehen und die DurchstoBpunkte durchsetzen 
zur Erzeugung wenigstens eines CT-Bildes Anordnung 
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